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I. INTRODUCTION 



Many-body states of nuclear shell model (SM) are linear combinations of Slater determi- 
nants built of bound single-particle (s.p.) states. Thus, the SM describes bound many-body 
systems which are isolated from the external environment of scattering states and decay 
channels. A natural generalization of the SM for weakly bound or unbound many-body 
systems, the so-called Gamow shell model (GSM), has been formulated recently in Berggren 
ensemble consisting of bound s.p. states, s.p. resonant states (Gamow [l| or Siegert ^ 
states), and the complex-energy s.p. continuum states 3|, |4|. A complete set of s.p. states 
can be defined in the Bergeren ensemble a a. The complete set of many-body states is 
then given by all Slater determinants spanned by s.p. states of the complete s.p. set in 
the Berggren ensemble 3, J]. So defined theoretical framework gives the full description of 
interplay between scattering states, resonances and bound states in the many-body wave 
function, without imposing a limit on the number of particles in the scattering continuum. 
On the other hand, the asymptotic decay channels in GSM are not individually resolved. 
Hence, this approach cannot be applied for a description of nuclear reactions and remains 
the tool for nuclear structure studies. 

The unification of structure and reaction theories is possible in the continuum shell model 
(CSM) formalism [ijsl 9L including the recently developed shell model embedded in the 
continuum (SMEC) [lO, U,Q, [3- In this formalism, the s.p. basis includes bound states 
and the real-energy continuum states. Feshbach's projection technique [jjj], used in CSM 
and SMEC, allows to describe on the same footing the nuclear reactions, including the 
rearrangement term, and the nuclear structure of well-bound, weakly-bound or unbound 
nuclear states. 

Bound and scattering s.p. states define two orthogonal subspaces go an d Po, respectively. 
The s.p. resonances do not belong to the Hilbert space and have to be regularized before 
including them in the SMEC framework. The regularization procedure consists of including 
s.p. resonances in a discrete part of the spectrum after removing the scattering tails which are 
included in the embedding scattering continuum. These regularized resonances are usually 
called the quasi-bound states embedded in the continuum (QBSEC). The new subspaces: 
q including bound and QBSEC s.p. states, and p including non-resonant scattering states 
and scattering tails of regularized resonances, are subsequently reorthogonalized. 
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The many-body states in SMEC are given by all Slater determinants spanned by s.p. 
states in q and p. Many-body states with all particles occupying bound and QBSEC s.p. 
states span the Q = <2 subspace of the Hilbert space. The complement subspace V includes 
many-body states with one or more particles in the scattering states, hence: V = Yli=i Qi> 
where Qj projects on the space with % particles in the continuum. 

The number of particles in the scattering continuum provides a natural hierarchy of ap- 
proximations in the CSM. Technical complications associated with inclusion of multiparticle 
continua and complex asymptotic channels are such that none of the early CSM or SMEC 
studies considered more than one particle in the continuum. Only recently, two-particle 
continuum has been included in the SMEC for the description of the two-proton radioactiv- 



ity 



131 ] . The general framework unifying the reaction theory and the structure theory for 



problems with any number of particles in the scattering continuum have been worked out 
as well {isl |. 

One of key elements in the practical implementation of SMEC scheme is the definition 
of orthogonal subspaces q, p and, consequently, Qo, Qi, Q2, ■ ■ ■ many-body subspaces. This 
definition is associated with the extraction of regularized resonances from the s.p. scattering 
continuum. In all previous SMEC applications (for an extensive list of applications see Refs. 
lol . H, 12, 13 1) as well as in CSM studies of Dresden group jl], 16], the Wang and Shakin 



method was employed [17]. In this method, the construction of a QBSEC depends on few 
parameters: the (real) energy of a scattering wave function which is then regularized and 
removed from the continuum, and the parameters of the cutting function (Heaviside or Fermi 
functions) chopping off the tail of this wave function. The so defined QBSECs are auxiliary, 
artificial objects which do not correspond to any solution of the Schrodinger equation. 

Unfortunately, the CSM/SMEC results concerning, e.g., the partial decay widths, depend 
on the parameters of the cutting function. This ambiguity cannot be totally removed even 



if the condition on the s.p. width is applied in choosing the cutting radius 



In practical 



applications, the radius of the cutting function is selected close to the top of the Coulomb 
, ll3[ what yields a sensible prescription for narrow s.p. resonances. Problems with 



barrier 



12 



171 ] do not have the correct 



the choice of cutting function and the extraction of QBSECs appear for broad resonances. 
Moreover, the QBSECs obtained using Wang and Shakin method 
asymptotic behavior. 

In this paper, we present the new method of regularizing s.p. resonances which is unam- 



3 



biguous, parameter-free and yields states with the correct bound-state asymptotics. These 
bound states embedded in the non-resonant continuum are called the anamneses of reso- 
nances in the space of square- integrable functions 9£ 2 -functions0. 

In Sect. HIl we provide a short discussion of different radial solutions of the Schrodinger 
equation which are used in the construction of a complete s.p. basis in SMEC. Qn unam- 
biguous determination of the resonance anamneses is presented in Sect. 1111} and Sect. IIVI 
is devoted to the discussion of consequences of the extraction of resonance anamneses on 
the non-resonant, real-energy continuum states. Together, bound s.p. states, anamneses of 
s.p. resonances and a regularized real-energy s.p. scattering states provide the complete s.p. 
basis in the Hilbert space. This basis can be used to obtain the complete many-body basis 
in CSM/SMEC studies. 

Applications of the new method are discussed in Sect. |VJ We shall present examples of 
the anamneses for resonances of different widths and angular momenta £, both for neutrons 
and protons. General features of the resonance anamneses, such as the energy dependence 
of both the root-mean-square (RMS) radius or the matching point between inner and outer 
solutions of the Schrodinger solution, will be analyzed as well. Finally, main conclusions will 
be given in Sect. IVII 

II. RADIAL SCHRODINGER EQUATION: THE GENERAL CONSIDERATIONS 

Let us consider a spherical potential V(r) describing qn interaction between a nucleon 
and a target nucleus. In the center of mass coordinates, the one-body radial wave function 
u(r) of the relative motion is the solution of the Schrodinger equation: 



where fi is the reduced mass and £ is the relative angular momentum. In the following, we 
shall omit the angular momentum index £ to simplify notations. 
We are interested in three kinds of solutions of Eq. (CQ) : 

• The scattering solutions, which form a continuum for real and positive momenta k. 
These solutions are regular at r = and asymptotically take a form: 




(1) 




(2) 
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where are irregular Coulomb wave functions or Hankel functions for protons and 
neutrons, respectively. C ± in ([2]) are some constants. 

• The bound state solutions, which compose a discrete set for imaginary and positive 
values k = k n . These solutions have an asymptotic of outgoing waves: 

u n (r) ~ Ak n rh + (k n r) . (3) 

• The resonance solutions, which constitute a discrete set for k = k T ^ s such that 
TZe{k T ^ s ) > 0, Xm{k r ^ s < 0) and TZe{k r ^ s ) > — Xm(k 7 n es ). These solutions have also an 
outgoing wave asymptotic. 

The asymptotic form of irregular Coulomb wave functions (Eqs. (J2]), ([3])) is: 

7 -=f(M-i) 

h±(kr) ~ e ±i(*T-„log(2fcr)) ,q 

kr 

where t] = fie 2 ZiZ 2 /h 2 k is the Sommerfeld parameter. 

The numerical integration for these three kinds of solution is carried out first by inte- 
grating the regular solution in the inner region [0, R], where R is some matching radius. On 
this interval, the solution of Eq. (CTJ is written as: 

u(r) = Cu r k e9 (r) , (5) 

where the regular solution u r k eg verifies: 

YiYa.r~ l ~ l Uu 9 {r) = 1 , (6) 

and C is a constant. Then, we integrate the Jost functions H k which are the solutions of 
Eq. (PQ) with the asymptotic form when r — > oo: 

H±{r) ~ krh^kr) . (7) 

In the outer region [R, oof, the scattering solutions are: 

u k {r) = C-H^r) + C+Ht{r) . (8) 

The continuity of the logarithmic derivative, which ensures the matching between the two 
solutions (Eqs. (jSJ) and (jHJ)), gives two conditions: 

Cu^{R) = C-H^{R) + C + H+{R) 

(9) 

c du k e9 {r) =c -dH k -(r) dH+M 



dr 



r=R dr 
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r=R dr 



r=R 



for three unknown constants and C. The third equation is provided by the normalization 
condition for scattering states: 

c+c- = JL . po) 

On the interval [R, oof, bound and resonance wave functions are written as: 

u(r) = Afl+(r) . (11) 

The matching between the internal ([3} and external (II ip solutions is ensured for those 
discrete values of k n which nullify Wronskian of the regular solution with the Jost solution 

Hp- 

W(u r k e9 ,H+) (r) = . (12) 

For resonance states, the values of k n correspond to the poles of the scattering matrix (S- 
matrix), whose matrix elements are given by the ratio of outgoing and incoming waves in 
the asymptotic expression <^}. 

Properties of the scattering continuum are very sensitive to the position of these poles in 
the complex /c-plane. For a narrow resonance, the scattering states with /c-values on the real 
axis just above the corresponding pole of the S- matrix are localized in the inner region of 
a potential and resemble bound states with, however, an oscillatory asymptotics instead of 
an exponential asymptotics expected for standard bound state wave functions. This generic 
feature of near-pole scattering wave functions helps to define the QBSEC function by cutting 



off the oscillatory tail with either the Heaviside function [17| or the Fermi [12j function. 



III. DETERMINATION OF THE ANAMNESES OF S.P. RESONANCES 

The significance of a QBSEC function in the CSM/SMEC theoretical framework is due 
to the similarity between a near-pole scattering wave function in the inner region (r < R) 
and a bound state wave function. The behavior of QBSEC in the outer region (r > R) 
does not correspond to any solution of the Schrodinger equation and, in this sense, can be 
considered unphysical. In the following, we shall show that one can define an alternative 
to QBSEC functions, the anamneses of s.p. resonances which belong to the space of C 2 - 
functions. These resonance anamneses are bound s.p. states in the scattering continuum. 
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They provide a natural continuation of weakly bound s.p. states for positive energies and 
extract all resonance features from the scattering continuum. 

Given a pole of the S'-matrix at k res in the complex fc-plane, this similarity is particularly 
striking for scattering states with k close to: 



k = y TZey{k res ) 2 j . (13) 

It is then quite natural to select the scattering state with a real eigenenergy e res = h 2 n 2 /2fi 
for a construction of the resonance anamnesis. 

In the inner region [0, R], where R is yet arbitrary matching radius for inner and outer 
solutions, the resonance anamnesis should be proportional to the regular solution of Eq. (CQ) 
with k = k: 

v K (r) = Cu^(r) , (14) 

where C is some constant. 

In the external region ([R, oof), we require that the anamnesis of the resonance has the 
bound state asymptotic corresponding to k — in. Using the Jost solution Hf K of Eq. (Tj[|) 
with k = in, we write the wave function of resonance anamnesis in the interval [R, oof as: 

v K (r)=AH+(r) , (15) 

where A is some constant. Since the wave number is fixed, both in the inner [0, R] and outer 
[R, oof regions, therefore the continuity of the resonance-anamnesis wave function and its 
first derivative at r = R is equivalent to the condition: 

WM?,H+)(r) =0 . (16) 

r=R 

Since the wave number is fixed, both in the inner [0, R] and outer [R, oof regions, Eq. ( fl6l) 
becomes a condition fixing the value of the matching radius r m = R. Hence, the wave 
function of the resonance anamnesis on [0, oof can be as: 



v K {r) = C 



< 9 {r) (l - 8(r - r m )) + |S^tf+(r)9(r - r 



- r m =R 

(17) 



in this equation denotes the Heaviside function and the constant C is fixed by the nor- 
malization condition for v K (r). 



The anamnesis of a s.p. resonance should be orthogonal to other bound and scattering 
states of a given r 2 , £, j. If a bound state exists with the same angular momentum £,j, then 
the wave function of the resonance anamnesis v K (r) (cf Eq. [T7I) has to be orthonormalized 
with respect to this state. Construction of an appropriate non-resonant scattering continuum 
will be discussed in Sect. HVl 

The above procedure yields for the resonance anamnesis a C(l)-function with an appro- 
priate bound state asymptotic. It is important to notice that the radius r m (= R) in (fTT|) 
is determined by solving Eq. (Tl6|) . and does not correspond to any arbitrary cutting radius 
as in the procedure of Wang and Shakin [^\ . We will see in Sect. |V]that despite the fact 
that r m can in principle take any value in the interval [0, oof, the above method provides 
the resonance anamneses which maintain localized aspects of the corresponding resonance 
wave functions. 



IV. CONSTRUCTION OF THE NON-RESONANT CONTINUUM 



The radial wave functions corresponding to resonance anamneses are not eigenf unctions 
of the original Hamiltonian h which is used to define s.p. bound and scattering states. 
It is however possible to redefine this Hamiltonian in such a way that bound states of h, 
resonance anamneses, and the states of a non-resonant scattering continuum are generated 
by one and the same Hamilton operator. 

Let us denote by {|w n )} the state vectors corresponding to bound states, and by {|f m )} the 
state vectors corresponding to the resonance anamneses defined in the previous section. The 
energy of each resonance anamnesis corresponds to the resonance energy: e^ s = h 2 K^ n /2^i, 
where K m is given by ffl3|) . and an index m enumerates different poles. Bound states and 
resonance anamneses form together a discrete subset {|tt n )} of the complete set of basis 



states in Hilbert space 18|, Il9j ]. 



Let us now define a new Hamilton operator h: 

h = ^ \u n )e n {u n \ + php , (18) 

n 

where p is a projection operator on the non-resonant s.p. continuum: 

P = 1 - \un)(u n \ , (19) 



and 

e n for a bound state 

(20) 

e^ es for an anamnesis of the resonance state 
All vectors in the set {|"S n )} are the eigenstates of h: 

(e n -h)\u n ) = . (21) 

The scattering continuum is renormalized by the second term on the r.h.s. of Eq. (fl8l) . 
Below, we shall discuss the construction of scattering states for the modified Hamiltonian h 
(cf Eq. ( ITS]) ). Such states, that we denote \u), have to be solutions of the equation: 

(e-php)\u) = . (22) 
Since \u) is orthogonal to all bound states including resonance anamneses, we can write: 

p(e-h)\u) = . (23) 
The general solution of Eq. (1221) can be written as: 

\u) = \u h )+J2 a n\Un) , (24) 
n 

where \u h ) is a solution of the homogeneous equation: 

(g-h)\u h ) = , (25) 

{|[/ n )} are solutions of the particular equations: 

(e-h)\U n ) = \u n ) , (26) 

and {a n } are constants. The orthogonality of \u) with respect to all bound states and 
anamneses of resonances: 

(u m \u) = , (27) 

can be expressed as: 

A + Ba = 6 , (28) 
where components of the vector A are given by: 

A m = (u m \u h ) , (29) 



and those of the matrix B by: 

B m n = (u m \U n ) . (30) 
Inverting the linear system of equations (|28|) . we obtain finally: 

\u) = \u h )-J2\Un)B^A m . (31) 

n,m 

The homogeneous solution has been discussed in Sect. [Til Its asymptotic behavior is given 
in Eq. ([2]). The states {|t/ n )} are solutions of inhomogeneous equations with outgoing wave 
asymptotics: 

U n (r)~C n krh+(kr) , (32) 

where {C n } are constants. Inserting (T5]) and fl32|) into the radial representation of Eq. fl3TT) . 
one derives the asymptotic behavior of the general solution: 

u{r) — C~ (_iy+ 1 e - i O' - -')log(2fcr)) _|_ £ e i(fcr-J?log(2fcr)) ^ ^33^ 

where the scattering matrix is given by: 

S = (C + ~ E C n B^ A m ) = , (34) 

n,m 

and 5 is the non-resonant phase shift. 

A possible test to see whether this procedure correctly removes the s.p. resonances from 
the s.p. scattering continuum is to compare the real-energy continuum phase shift for both 
Hamiltonians: h and h. One demands that resonances are suppressed from the scattering 
continuum of h. On the other hand, the background phase shift calculated for h and h 
should be essentially the same. 



V. DISCUSSION OF THE RESULTS 



In this section, we shall discuss examples of resonance anamneses constructed by the 
method described in Sects. [TTJ and IIHI We shall also analyze properties of an associated 
non- resonant scattering continuum described in Sect. IIVI 

The starting point is the construction of an initial s.p. basis containing bound and 
scattering states. This basis is generated by the spherical Woods-Saxon (WS) potential 
consisting of both central and spin-orbit parts, and a Coulomb potential: 

V(r) = -V f(r) - 4V/ so (/ • s) 1 -^ + V c {r) , (35) 
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where 



f(r) 



1 + exp 



Rr 



d 



-i -1 



(36) 



In all examples, the WS potential has the radius i?o=3.5 fm, the diffuseness d =0.5 fm, 
and the spin-orbit strength V^ =3.5 MeV. The mass of a target nucleus is 16 amu. The 
Coulomb potential V c is assumed to be generated by a uniformly charged sphere of radius 
Ro and charge number Z=8. The depth of the central part is varied to simulate different 
situations. The resonance anamneses are extracted from the scattering continuum of the 
WS potential (135]) . Bound states, anamneses of resonances and remaining continuum states 
are subsequently reorthogonalized. 



A. Radial wave functions of resonance anamneses and associated phase shifts 

As an illustration of the method, we shall consider the 0d 5 / 2 proton resonance in the 
WS potential fl35|) with Vo=40 MeV. The 0d 5 / 2 resonance pole is situated at TZe(k res ) = 
0.265 fm -1 , Tm(k res ) = — 0.0014 fm -1 , what corresponds to the energy e res =1.55 MeV and 
the width r res =32.692 keV. The radial wave function of a scattering state with real energy 




10 20 30 40 50 
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FIG. 1: The radial wave function of a 0(5/2 proton scattering state at an energy corresponding to 
the real part of a 0(f5/2 resonance state energy. For more details, see the description in the text. 



e res is presented in Fig. [TJ One can see the localization of the wave function in the inner 
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region (up to ~10 fm). Solving the matching condition (fT6l) yields the matching radius: 
r m (= R)=18 fm. 

0,6 




"0 
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FIG. 2: The resonance anamnesis corresponding to the proton scattering state shown in Fig. 

m 



The resonance anamnesis corresponding to the proton resonance (see Fig. [2]) is 

extracted from the d§/2 scattering wave function function at the energy e res corresponding 
to the real part of the resonance energy (cf Fig. [T]). One may notice a hump at 

larger distances (10 fm < r < 20 fm) which is a characteristic feature of many resonance 
anamneses. 

^5/2 phase shifts of the resonant continuum generated by the Hamiltonian /iws with the 
WS potential fl35l) . and of the non-resonant continuum generated by the modified Hamilto- 
nian h ws (cf Sect. HV|) . are plotted in Fig. [3J One can see that the resonant part of phase 
shift is fully removed by the G?5/2-resonance anamnesis and the phase shift calculated for the 
scattering continuum of /iws exhibits only a smooth energy dependence. 

Another examples of the radial wave functions of resonance anamneses and the respective 
phase shifts are shown in Figs. HH3 Fig. H] shows the resonance- anamnesis wave functions 
corresponding to lpi/2 neutron resonances at three different energies. The depth of the 
WS potential is varied to yield the resonance energy e res = 0.036 MeV (the bottom part), 
0.112 MeV (the middle part), and 0.47 MeV (the upper part). The widths of these resonances 
vary from ~ 10 keV to ~ 560 keV. The resonance anamneses are extracted fromp!/2 neutron 
scattering wave functions at energies e res corresponding to the real part of the lpi/2 neutron 
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FIG. 3: Resonant (full line) and non-resonant (dashed line) d^j2 phase shift as a function of the 
energy e = h 2 k 2 /2fi in the center of mass. 



resonance energy. In all these three cases, well pronounced hump at large distances can be 
seen. 

Phase shifts in resonant and non-resonant scattering continua corresponding to those 
three lpi/2 neutron resonances shown are plotted in Fig. [5j The resonant continuum is 
generated by the /iws Hamiltonian, whereas the non-resonant continuum is given by the 
modified Hamiltonian /iws- O ne can see that the regularization procedure fully subtracts 
the resonant part of the phase shift, both for broad and narrow resonances, and the phase 
shift corresponding to the continuum of /iws exhibits only a smooth energy dependence. 

Radial wave functions of resonance anamneses and corresponding phase shifts for lsi/2 
proton resonance at three different energies are shown in Figs. |6] and [7J The depth of the 
WS potential is varied to yield the resonance energy at e res = 0.53 MeV (the bottom part), 
1.11 MeV (the middle part), and 1.29 MeV (the upper part). The resonance width in these 
cases vary from ~ 20 keV to ~ 600 keV. Also in this case, the extraction of the resonance 
anamnesis removes all resonant features from the scattering continuum of hws- 

In all cases, the anamnesis of a resonance constructed by the method proposed in this 
work provides an excellent description of localized aspects of the resonant wave function 
in a broad range of resonance widths. The so defined resonance anamnesis becomes an 
image of the resonance wave function in the subspace of £ 2 -functions. Its removal from the 
subspace of scattering wave functions leaves only the non-resonant continuum with smoothly 
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FIG. 4: Radial wave functions of resonance anamneses corresponding to P1/2 neutron scattering 
states at three different energies of lpi/2 neutron resonance. For more details, see the description 
in the text. 

varying phase shifts. This feature of resonance-anamnesis wave functions give them a special 
significance in a decomposition of the scattering continuum. 

B. RMS radius and matching radius of the resonance-anamnesis wave functions 

In the following, we shall investigate the energy dependence of the resonance-anamnesis 
radial wave functions by calculating the RMS radius and the matching radius. 



14 





2 r 




1- 












0- 








-1 ■ 




-4 





p res _ n A 7 A/T P V 




r = 557.8 keV 


\ 

\ 

\ 

\ 









-3> 

2 



" 2 o L 



0,5 1 1,5 2 

e (MeV) 



e res = 0.112 MeV 
T = 58 keV 



0,2 0,4 

e (MeV) 



e res = 0.036 MeV 
T = 9.96 keV 



0,1 0,2 

e (MeV) 



2,5 



0,6 



FIG. 5: Resonant (full line) and non-resonant (dashed line) p\/2 neutron phase shift for three 
resonances shown in Fig. [H 



The dependence of the mean-square radius of a bound state on the distance e to the 



continuum threshold has been studied by Riisager et al 
et al. 



201 ] for spherical nuclei and Misu 



2l| for deformed nuclei. (For the discussion of three-body halo asymptotics, see Ref. 
22j.) The mean-square radius of a spherical neutron orbit varies as (— e) _1 for I = and 
as (— e)" 1 / 2 for t — 1, but remains finite for higher angular momenta due to the centrifugal 
barrier that confines the radial wave function. For proton bound states, the mean-square 
radius is finite for all i at the continuum threshold because of the Coulomb barrier. 

The resonance-anamnesis wave functions constructed by the method proposed in this 
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FIG. 6: Radial wave functions of resonance anamneses corresponding to S1/2 proton scattering 
states at three different energies of ls]y 2 proton resonance. For more details, see the description 
in the text. 



work have an asymptotic behavior of a bound state wave function with momentum 
■J TZe(^(k res ) 2 ^j . Consequently, it is likely that certain near-threshold features of resonance- 
anamnesis wave functions are similar to features of weakly bound s.p. states. In other 
words, one expects not only that the above method provides states which provide a faithful 
image of localized aspects of resonances, but also that it provides states which are a natural 
continuation of weakly bound s.p. states into the continuum. 

Figs. M and [9] present the RMS radius of a weakly bound s.p. state and a corresponding 
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FIG. 7: Resonant (full line) and non-resonant (dashed line) s\ji proton phase shift for three 
resonances shown in Fig. [6J 



resonance anamnesis as a function of the energy for lpi/2 neutron and lsi/2 proton s.p. 
orbits, respectively. The curve for negative energies corresponds to the RMS radius for 
bound s.p. states: neutron lpi/2 (Fig. [H]) and proton lsi/ 2 (Fig. [§]) states. For positive 
energies, the curve shows the RMS radius of the corresponding resonance-anamnesis wave 
function extracted from the scattering wave function at the energy which is equal to the 
real part of the resonance energy, respectively. The corresponding bound and real-energy 
scattering states are found by varying the depth of the WS potential fl35|) . For neutron lpi/2 
orbits (cf Fig. [8]), the depth of the central part Vq varies from 69.012 MeV to 61.44 MeV, 
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FIG. 8: RMS radius of lpi/2 neutron orbit as a function of the energy. The branch for positive 
energies is obtained for lp\/2 resonance anamnesis constructed by the method described in Sect. 
IIIII For more details, see the description in the text. 
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FIG. 9: RMS radius of the lsi/2 neutron orbit as a function of the energy. The curve for positive 
energies is obtained for the lsi/ 2 resonance anamnesis constructed by the method described in 
Sect. IIIII For more details, see the description in the text. 
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whereas for proton lsi/2 orbits, Vq changes from 50.233 MeV to 39.06 MeV. 

The gap seen in Figs. [8] and [9] close to e = comes essentially from numerical limitations 
in the construction of reliable wave functions for a weakly-bound state and its resonance 
anamnesis at positive energies. One may notice a slight asymmetry in the energy dependence 
of RMS radius on both sides of e = 0. This is essentially an effect of the centrifugal and/or 
Coulomb barrier for states with e > 0. 

In all previous CSM or SMEC studies [lO, 



11 



12 



131]). the QBSEC radial wave function 



was taken to be proportional to a regular solution of Eq. (TTJ) up to a certain radius r cut , 
and for r > r cut this solution was cut by a cutting function. The value of r cut in such an 
approach is a free parameter. For most of CSM/SMEC applications, the value of r cut was 
close to the position of the top of the centrifugal and/or Coulomb barrier. 
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FIG. 10: Matching radius r m (full line) in function of the energy e of the resonant lpi/2 neutron 
state. The depth Vo of the potential is varied from 66.1 MeV to 61.44 MeV. For comparison, the 
external turing point (dashed line) is presented too. 



As discussed in Sect. the resonance-anamnesis wave functions which replace QBSEC 
functions, are proportional to the regular solution of Eq. (CQ) up to the matching radius 
r m = R. For r > r m , the resonance anamnesis takes an asymptotic form of a bound state 
with a well defined wave number. The matching radius is not a free parameter of the theory 
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but its value is uniquely determined by solving Eq. (fT6|) for r. At r = R, inner and outer 
solutions join smoothly to form a C(l) radial wave function of the resonance anamnesis. 

The matching radius r m , whose position depends on features of the studied resonance, is 
an interesting supplementary information about the resonance-anamnesis wave function. r m 
depends uniquely on resonance characteristics, such as the energy (width), angular momen- 
tum and t z . The energy (width) dependence of r m is generic, depending only on the angular 
momentum and r z . This dependence is illustrated in Fig. [10] on an example of neutron lpi/2 
wave function of the resonance anamnesis. We can see that r m grows with decreasing e. 
The dashed line in Fig. [TO] shows also the variation of an external turning point r etp of the 
centrifugal barrier with the resonance energy. Except for a shift, both functions r m (e) and 
r e*p(e) are rather similar, in particular for large e. In all studied cases, the matching radius 
lies outside of the external turning point of the potential barrier. 



VI. CONCLUSIONS 



Construction of the many-body basis in Hilbert space is a key problem in the SM de- 
scription of open quantum systems, such as the weakly bound nuclei, radioactive decays, or 
low-energy nuclear reactions. This construction is based on the determination of a complete 
s.p. basis which consists of discrete states (subspace q) and a scattering continuum (subspace 
p). A consistent unified formulation of nuclear structure and reactions is possible within the 
CSM/SMEC if the s.p. resonances are removed from the scattering continuum p to be put 
in q [8j. This regularization procedure for s.p. resonances is associated with an extraction 
of localized part of s.p. resonances from p which leaves only non-resonant scattering states. 

In this paper, we have developed an unambiguous, parameter-free method of extracting 
the localized component of the s.p. resonance which allows to determine a new subspace 
q of discrete s.p. states, consisting of bound states and resonance anamneses, and a new 
subspace p of non-resonant scattering states. The resulting s.p. basis is complete and can 
be used for the construction of complete many-body basis in the Hilbert space. This new 
s.p. basis could also provide an interesting alternative for solving Hartree-Fock-Bogolyubov 
equations with finite-range density-dependent interactions for weakly-bound systems. 

The resonance anamneses provide a faithful image of localized features of resonances in 
the space of £ 2 -functions, removing all resonant aspects from calculated phase shifts. These 
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anamneses of the resonance states provide also a smooth continuation of weakly-bound s.p. 
states into the low-energy continuum. Appealing features of those states give a special 
significance to the s.p. basis which includes them in its discrete part. Future applications 
of this particular basis may open new and yet unexplored horizons for a direct calculation 
of multi-particle resonant states in nuclear systems. It may also be useful in extracting a 
resonance contribution from many-body observables in weakly-bound or unbound systems. 
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